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Double perovskite Bi2FeCr06, related with BiFeOs, is very interesting because strong ferroelec- 
tricity and high magnetic Curie temperature beyond room temperature are observed in it. However, 
existing density-functional-theory (DFT) studies, using pseudo-potentials, produce metallic ground 
state under the local density approximation (LDA) and need LDA+U method to yield needed non- 
metallic ground state, resulting in low magnetic Curie temperature (below 130 K). Here, we optimize 
its crystal structure and then investigate its electronic structure and magnetic and optical proper- 
ties by combining the full-potential augmented plane wave method with Monte Carlo simulation. 
Our optimized structure is a robust ferrimagnetic semiconductor. This nonmetallic phase is formed 
due to crystal field splitting and spin exchange splitting, in contrast to Mott-Hubbard states in 
previous DFT studies. Spin exchange constants and optical properties are calculated. Our ab initio 
magnetic Curie temperature is 450 K, much higher than previous DFT-based value and consistent 
with experimental results. Our study and analysis reveals that the main magnetic mechanism is an 
antiferromagnetic superexchange between Fe and Cr over the intermediate O atom. These results 
are useful to understanding such perovskite materials and exploring their potential applications. 

PACS numbers: 75.50.-y, 71.20.-b, 75.10.-b, 85.75.-d 



I. INTRODUCTION 

Magnetic materials keeping high spin polarization at 
room temperature or higher are highly desirable for spin- 
tronic applications^ 2 -. Half-metallic materials 3 ' and fer- 
romagnetic (or ferrimagnetic) semiconductors are inten- 
sively studied for this purpose. Up to 96% spin polar- 
ization has been achieved experimentally in the cases of 
(La,Sr)Mn03 and Cr02 materials 2,4 . Since 1998, double 
perovskite oxides have been extensively explored. Both 
half-metals and semiconductors with macroscopic mag- 
netization beyond room temperature have been found in 
such materials 5 9 . High quality samples of semiconductor 
technology standard have been achieved for some of these 
compounds 10 . Furthermore, it has been established that 
ferroelectricity and magnetism can coexist above room 
temperature in BiFeOs materials^"— . These stimulate 
world-wide interest in exploring multiferroic materials for 
novel devices^"—. 

Double perovskite Bi2FeCr06 can be constructed from 
BiFeOs, and is first predicted to have both ferroelectric- 
ity and ferrimagnetism, with theoretical electric polariza- 
tion 80 /iC/cm 2 and magnetic Curie temperature below 
130 Ki£i±£. Unfortunately, the Mott-Hubbard ground 
state excludes higher Curie temperature in such theo- 
retical studies. Subsequently, the Bi2FeCr06 is synthe- 
sized successfully by several groups, and experimental 
measurements on good samples reveal electric polariza- 
tion up to 60 /iC/cm 2 and magnetic Curie temperature 
well above room temper at ureii - —. The best experimen- 
tal Curie temperatures are much higher than the first- 
principles predicted one, while the electric polarization 
is consistent between theory and experiment. Reason- 



able explanation of this contradiction is highly desirable 
to understanding the basic physics of such materials and 
exploring their applications. 

In this paper, we combine a full-potential density- 
functional-theory (DFT) method with Monte Carlo sim- 
ulation to investigate double perovskite Bi2FeCr06, with 
emphasis on its electronic structure and magnetic Curie 
temperature. We shall first optimize its structure with 
GGA and then study its electronic structure. We 
shall show that the resultant ferrimagnetic semiconduc- 
tor phase is robust against both change of exchange- 
correlation potential and variation of structural parame- 
ters. Furthermore, we shall calculate spin exchange con- 
stants and thereby determine the magnetic Curie tem- 
perature. Our ab initio magnetic Curie temperature is 
450 K, much higher than previous DFT-based value but 
consistent with experimental results. In addition, we 
shall present calculated optical properties. The magnetic 
mechanism is understood in terms of the electronic prop- 
erties. More detailed results will be presented in the fol- 
lowing. 

The rest of this paper will be organized as follows. In 
the next section, we shall describe our computational de- 
tails. In the third section we shall describe first-principles 
structural optimization and present our optimized re- 
sults. In the fourth section we shall present main first- 
principles electronic structures. In the fifth section we 
shall present first-principles magnetic calculations and re- 
sulting Curie temperature. In the sixth section we shall 
present first-principles optical properties. In the seventh 
section we shall make some necessary discussions. Fi- 
nally, we shall give our conclusion in the eighth section. 



II. COMPUTATIONAL DETAILS 

We use the full-potential augmented plane wave 
method within the density functional theory (DFT) 24 , 
as implemented in package WIEN2k 25 . Generalized gra- 
dient approximation (GGA) 26 is used for the exchange- 
correlation potential for calculating our main results, 
and local density approximation (LDA) 27 and a modi- 
fied Becke- Johnson exchange potential (mBJ) 28 are also 
used for comparative calculations. Full relativist ic effects 
are calculated with the Dirac equations for core states, 
and the scalar relativistic approximation is used for the 
other states 2 ^"—. The spin-orbit coupling is neglected 
because it has little effect on our main conclusion (to 
be detailed in the following). The muffin-tin radii of 
Bi, Fe, Cr, and O are set to 2.13, 1.84, 1.84, and 1.64 
bohr. The spherical expansion in the muffin-tin spheres 
is done up to l max =10. The parameter R mt x K max is 
set to 7.5. We use 1500 k-points in the first Brillouin 
zone (226 k-points in the irreducible Brillouin zone) dur- 
ing the calculations. The self-consistent calculations are 
considered to be converged only when the integration of 
absolute charge-density difference per formula unit be- 
tween the successive loops is less than 0.0001 |e|, where e 
is the electron charge. 

Metropolis algorithm and variants are used for our 
Monte Carlo simulations^ 2 -^. Several three-dimensional 
lattices of up to 30x30x30 NaCl unit cells with periodic 
boundary condition are used in these calculations. The 
first 50,000 Monte Carlo steps (MCS) of total 100,000 
MCS are used for the thermal equilibrium, and the re- 
maining 50,000 MCS are used for a given temperature. 
The T c value is determined through investigating the av- 
erage magnetization and magnetic susceptibility as func- 
tions of temperature^. 



III. STRUCTURAL OPTIMIZATION 

It has already been prove d 15 i 16 that the lowest energy 
structure of double perovskite Bi2FeCr06 has R3 symme- 
try (space group # 146). Its structural difference from 
the cubic double perovskite is mainly due to the alternat- 
ing rotations of oxygen octahedra around Fe/Cr cations. 
The magnetic order turns out to be ferrimagnetic with 
every Cr spin orienting upward and every Fe spin down- 
ward. We use experimental lattice constant a = 5.537 A 
and c = 13.502 A (a = 60.20°) which are intermediate 
values of existing experimental results 17 23 . After inter- 
nal structure optimization, every atom has a displace- 
ment from that of the cubic perovskite structure. We 
show its crystal structure in Fig. 1. The optimized pa- 
rameters are summarized in Table HI The bond lengthes 
and Fe-O-Cr bond angle can be calculated from the struc- 
tural parameters. The calculated values are summarized 
in Table HI 

These structural parameters are not much different 
from the previous LDA-optimized result s 15 1 16 . Actually, 
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FIG. 1: (color online) The Rc (#146) Crystal structure of 
double perovskite Bi2FeCr06- The right-side is the primitive 
cell, and the left-side is used to describe the relative positions 
of the atoms. O, Fe, Cr, and Bi are denoted by colorful balls 
with different size from the smallest to the biggest. 



TABLE I: Optimized internal structural parameters of double 
perovskite Bi2FeCr06 with the Rc (#146) crystal structure. 



Atom 


Parameter 


Values 


Bi 


Zl 


0.0005, 0.5040 


Fe/Cr 


Z2 


0.7330/0.2260 


O 


X 


0.5433, 0.0511 




y 


0.9542, 0.9043 




z 


0.3937, 0.4459 



we obtain 55.2/iC/cm 2 for the sum s ^ i ZiUi of formal 
charge Z{ times displacement using formal charges 
3.0, 3.0, 3.0, and -2.0 for Bi, Fe, Cr, and O, respectively. 
This value is almost the same as the previous 55/iC/cm 2 
in Refs. 15 and 16. Therefore, we believe that simi- 
lar ferroelectric polarization can be calculated with the 
present structural parameters. The total magnetic mo- 
ment is equivalent to 2/i# per formula unit. Different 
magnetic moments, 3.70/i# and 2.18/i^, are observed in 
the spheres of Fe and Cr atoms, and magnetic moments 
in the spheres of other atoms are tiny. These implies that 
Fe and Cr are both in the high spin states. 



TABLE II: Bond lengths (A) and Fe-O-Cr angle (°) of double 
perovskite Bi2FeCr06 with the Rc (#146) crystal structure. 



Bond length 


Fe-0 


2.130, 1.971 




Cr-0 


1.947, 1.994 




Bi-0 


2.347, 2.387 






2.406, 2.416 


Bond angle 


Fe-O-Cr 


153.7, 153.6 
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IV. ELECTRONIC STRUCTURE 

With the optimized crystal structure, we can study 
the electronic structure of the E^FeCrOe- In Fig. 2 we 
present the spin-resolved total density of states (DOS) 
and partial DOS projected in the d states of Fe and Cr be- 
tween -7.3 and 5 eV. The corresponding spin-dependent 
band structure is presented in Fig. 3. It is clear that 
the Bi2FeCr06 is a ferrimagnetic semiconductor, in con- 
trast with the ferrimagnetic metallic phase obtained in 
previous pseudo-potential calculation o 15 i 16 i 34 i 35 . 




Energy (eV) 

FIG. 2: (color online) Spin-resolved density of states (DOS, in 
state/eV per formula unit) of double perovskite Bi2FeCr06, 
calculated with GGA. The solid line is total DOS, and dashed 
and dotted lines refer to partial DOS projected in the atomic 
spheres of Fe and Cr, respectively. The upper part in each 
panel is majority-spin DOS result, and the lower the minority- 
spin one. 




FIG. 3: (color online) Spin-resolved energy bands of double 
perovskite Bi2FeCr06, calculated with GGA. The left-side 
panel is for majority-spin, and the right-side for minority- 
spin. 

In the majority-spin channel, the bands between -7.3 
and -1.4 eV are from the 18 O 2p states; the three bands 
between -0.8 and eV originate mainly from Cr 3d t2g; 
the three empty bands are almost Fe 3d t2g states; and 
the upper states are from Fe 3d eg, Cr 3d eg, and Bi 



6p. In the minority-spin channel, the 23 filled bands are 
from 18 O 2p states and 5 Fe 3d ones; and the empty 
bands from Cr 3d and Bi 6p states. Here, the Fe 3d 
states hybridize with the O 2p and the Cr 3d ones interact 
strongly with the Bi 6p, in contrast to the majority-spin 
channel. 

Because it is well known that both LDA and GGA un- 
derestimate semiconductor gap, we use a new exchange 
potential, mBJ, to improve the electronic structure, es- 
pecially the semiconductor gap. The mBJ DOS are pre- 
sented in Fig. 4. The main features are similar to the 
GGA results, except that the mBJ gap is much wider 
than the LDA and GGA ones. The three semiconductor 
gaps are summarized in Table HTT1 
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FIG. 4: (color online) Spin-resolved density of states (DOS, in 
state/eV per formula unit) of double perovskite Bi2FeCr06, 
calculated with mBJ. The solid line is total DOS, and dashed 
and dotted lines refer to partial DOS projected in the atomic 
spheres of Fe and Cr, respectively. The upper part in each 
panel is majority-spin DOS result, and the lower the minority- 
spin one. 



TABLE III: The semiconductor energy gaps (in eV) of double 
perovskite Bi2FeCr06, calculated with our optimized struc- 
ture and previous on o 15 ' 16 in terms of LDA, GGA, and mBJ. 



Structure 


LDA 


GGA 


mBJ 


this 


0.26 


0.42 


1.48 


previouo 15lls 


0.30 


0.48 





Our study shows that this feature of ferrimagnetic 
semiconductor is robust to a changing of crystal struc- 
ture parameters. For example, the GGA and LDA gaps 
in the lower line of Table IIIII are calculated with the pre- 
vious structure parameter o 15 i 16 . 



V. MAGNETIC CURIE TEMPERATURE 

In order to study spin exchange interactions, we com- 
pare the total energies of the four magnetic configurations 
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in the supercell consisting of 20 atoms, namely ferrimag- 
netic (the ground state) and ferromagnetic orders, and 
two other magnetic orders constructed by reversing one 
of the two Cr or Fe spins in the ferrimagnetic order. With 
respect to the ground state, the three magnetic structures 
have total energies: 310, 161, and 166 meV per formula 
unit, respectively. Because the induced spin density at 
the other atoms are very small compared to those at the 
magnetic atoms, we consider only Fe and Cr. The mag- 
netic energies reflect the spin exchange energies e^- where 
i and j denote two different spins. Fe and Cr atoms form 
a lattice of the NaCl crystal structure5*L2*JL2&. While the 
magnetic moments in the spheres of Fe and Cr are 3.70 
and 2.18 the Fe 3+ and Cr 3+ cations contribute 5 and 
3 hb, respectively. We can assign spin values s = § and 
s = | to the Fe and Cr spins, respectively. Therefore, we 
obtain the effective spin Hamiltonian: 



H 



^1 Jij &i 

(ij) 



(1) 



where Si is spin operator at site i (in both of the Fe and 
Cr sublattices), the summation is over spin pairs, and 
the spin exchange constant Jij is limited to the near- 
est (Fe-Cr) and the next nearest (Fe-Fe and Cr-Cr) spin 
pairs. There exists a relation between the energies e^- 



and constants J^, = JijSiSj, where Si takes either 2 



or |. The spin exchange energies and constants 
are summarized in Table IIVI 



TABLE IV: The calculated spin exchange energies (e^) and 
constants (Jij) of double perovskite Bi2FeCr06 for the nearest 
and next nearest spin pairs ((ij)), and the resulting Curie 
temperature T c from Monte Carlo simulation. 



(ij) 


Fe-Cr 


Fe-Fe 


Cr-Cr 


dj (meV) 


25.8 


-1.88 


-1.06 


Jij (meV) 


6.89 


-0.30 


-0.47 


Tc (K) 


450 







We carry out Monte Carlo simulations of the Heisen- 
berg model to estimate the T c of the materials 32,33 . It is 
well known that Curie temperature will be a little under- 
estimated if classical approximation to the Heisenberg 
model is used in the Monte Carlo simulation. We use 
Monte Carlo simulation of the classical Heisenberg model 
to give a lower bound for the Curie temperature of the 
quantum spin model (1). Our simulated average mag- 
netization (in unit of 2/i# per Fe-Cr pair) as a function 
of temperature from classical Heisenberg model is pre- 
sented in Fig. 5. Because the classical approximation 
is used, the low-temperature part is not very accurate 
with respect to the model (1), but the high-temperature 
part should be accurate enough to get reasonable Curie 
temperature. The calculated T c value is equivalent to 
450 K, meaning that real Curie temperature should be 
a little higher than 450 K. This Curie temperature is 
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FIG. 5: (color online) Temperature dependent magnetiza- 
tion (in unit of 2/jlb per Fe-Cr pair) of double perovskite 
Bi2FeCr06, simulated with Monte Carlo method and the clas- 
sical Heisenberg model. 



much higher than the original prediction (lower than 130 
K) 15 i 16 , but it is consistent with experimental high Curie 
temperatures beyond room temperatur o 17 i 2Q ~ — . 



VI. OPTICAL PROPERTIES 

Because optical properties are important to semicon- 
ductor, we calculate the direct electronic contributions of 
the optical dielectric function e(u) and optical conductiv- 
ity a(u) curves of the Bi 2 FeCr06 as functions of photon 
energy (uj). Considering that the semiconductor gaps 
have important impact on the low-energy parts of these 
optical functions, we calculate them with both GGA and 
mBJ. These results, both real (Re) and imaginary parts 
(Im) between and 12 (or 13) eV, are presented in Figs. 
6 and 7. 
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FIG. 6: (color online) Optical dielectric functions, both real 
(solid) and imaginary (dotted) parts, of double perovskite 
Bi 2 FeCr0 6 , calculated with GGA (thick) and mBJ (thin). 

For Re(e), the GGA result of the low energy limit is 
substantially larger than the mBJ one, both have simi- 
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FIG. 7: (color online) Optical electric conductivity, both real 
(solid) and imaginary (dotted) parts, of double perovskite 
Bi 2 FeCr0 6 , calculated with GGA (thick) and mBJ (thin). 
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FIG. 8: (color online) The majority-spin part of the Bi 6s, O 
2s, and Bi 5d energy bands of double perovskite Bi2FeCr06 
in the R3 (the left part) and cubic (the right part) structures, 
calculated with GGA. 



lar trend when the photon energy is larger than 1.5 eV. 
For Im(e), main difference appears for low energy (0-1.5 
eV). There is a sharp peak at 0.75 eV for the GGA re- 
sult, which is related to the steep decrease in the GGA 
curve of Re(e) at the same photon energy. The mBJ 
result of Im(e) is equivalent to zero in this region, reflect- 
ing the broad semiconductor gap from mBJ. For higher 
energy, both of the results are similar. For the ct(uj) 
curves, the main difference also appears in the low en- 
ergy region, Re(a(u)) having a peak approximately at 
1.0 eV and Im(cr(u;)) having a valley approximately at 
0.5 eV. For the higher energy, both of the <j(uj) curves 
have similar behaviors. These calculated optical func- 
tions could be useful in exploring the optical properties 
of the Bi2FeCr06 materials and others similar. 



VII. DISCUSSIONS 

Our DFT investigation with GGA, LDA, and mBJ 
shows that the ground-state phase of double perovskite 
Bi2FeCr06 is a ferrimagnetic semiconductor with a clear 
semiconductor gap (larger than 0.26 eV), in contrast to 
previous ferrimagnetic metallic phases from LDA pseudo- 
potential calculation ^ 15 ! 16 ! 34 ^ . Even using the previous 
structural paramter o 15 i 16 , we still obtain the ferrimag- 
netic semiconductor phase with both GGA and LDA. 
Furthermore, our spin exchange energies and interaction 
constants produce the high Curie temperature 450 K 
which is consistent with experimental high Curie temper- 
ature beyond room temperatur e - 17 ! 20- — . This is in con- 
trast to the previous LDA+U calculations which yield 
Mott- Hubbard insulator phases and low Curie temper- 
atures (below 130 K) 15 i 16 . These imply that we need 
not use LDA+U or GGA+U to achieve the reasonable 
ferrimagnetic semiconductor phase for double perovskite 
Bi 2 FeCr0 6 . 

Here, our semiconductor gap is formed between the 



filled Cr 3d t2g and the empty Fe 3d t2g bands. Es- 
sentially, it is originated mainly from the crystal field 
splitting due to the deformation of the O octahedrons 
plus the spin exchange splitting of the 3d electrons. Our 
calculated spin exchange interaction between the nearest 
Fe and Cr atoms is positive, and the next nearest (Fe-Fe 
and Cr-Cr) ones are negative. These spin interactions 
become substantially weaker from the nearest to next 
nearest neighboring spin pairs, which implies that such 
a spin description is reliable. The main spin interaction 
is intermediated by the O atom in between the nearest 
Fe and Cr atoms. The Fe and Cr spins contribute op- 
posite magnetic moments so that the ferrimagnetism is 
formed. Therefore, the main magnetic mechanism is the 
antiferromagnetic superexchange over the O atom. 

Now we address possible contributions from deeper en- 
ergy levels. For this purpose, we present in Fig. 8 the 
bands of Bi 6s, O 2s, and Bi 3d states of the R3 (space 
group #146) and cubic (space group #225) structures. 
They are below -9 eV from the Fermi levels. The Bi 5d 
bands (-23.3 or -24.4 eV) in both of the structures are 
completely flat, meaning that Bi 5d states are isolated, 
but there is a clear correlation between the two Bi 6s 
bands (from -11.0 to -9.4 eV) and two of the six O 2s 
ones (from -19.9 to -17.7 eV) in the case of the R3 struc- 
ture. Noticeably, its Bi 6s bands are more separated from 
the O 2s ones than those of the cubic structure. This 
can be attributed to the expelling due to the small but 
nonzero hybridization of the two sets of s states. The 
interaction between the Bi 6s and O 2s states should be 
involved in the ferroelectric property in the R3 structure, 
because they are absent in the cubic structure without 
ferroelectricity. 



VIII. CONCLUSION 

In summary, we have investigated the electronic struc- 
ture and magnetic and optical properties of double per- 
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ovskite Bi 2 FeCr06 by combining the full-potential aug- 
mented plane wave method with Monte Carlo simula- 
tion. Our calculations show that our full-potential GGA- 
optimized phase is a ferrimagnetic semiconductor and it 
is robust against both change of exchange-correlation po- 
tential and variation of structural parameters. This non- 
metallic phase is formed due to crystal field splitting and 
spin exchange splitting, in contrast to Mott-Hubbard in- 
sulating state in previous DFT studies. Spin exchange 
constants and optical properties are calculated. Our ab 
initio magnetic Curie temperature is 450 K, much higher 
than previous DFT-based value but consistent with ex- 
perimental results. Our study and analysis reveals that 



the main magnetic mechanism is an antiferromagnetic 
superexchange between Fe and Cr over the intermediate 
O atom. Our results are useful to understanding such 
materials and exploring their potential applications. 
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